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We investigate the thermal photon production-rates using one dimensional boost-invariant second 
order relativistic hydrodynamics to find proper time evolution of the energy density and the temper- 
ature. The effect of bulk- viscosity and non-ideal equation of state are taken into account in a manner 
consistent with recent lattice QCD estimates. It is shown that the non-ideal gas equation of state 
i.e e — 3P 7^ behaviour of the expanding plasma, which is important near the phase-transition 
point, can significantly slow down the hydrodynamic expansion and thereby increase the photon 
production-rates. Inclusion of the bulk viscosity may also have similar effect on the hydrodynamic 
evolution. However the effect of bulk viscosity is shown to be significantly lower than the non-ideal 
gas equation of state. We also analyze the interesting phenomenon of bulk viscosity induced cav- 
itation making the hydrodynamical description invalid. We include the viscous corrections to the 
distribution functions while calculating the photon spectra. It is shown that ignoring the cavitation 
phenomenon can lead to erroneous estimation of the photon flux. 



I. INTRODUCTION 

Thermal photons emitted from the hot fireball created 
in relativistic heavy-ion collisions is a promising tool for 
providing a signature of quark-gluon plasma [T]-[6] (see 
[343 for recent reviews). Since they participate only 
in electromagnetic interactions, they have a larger mean 
free path compared to the transverse size of the hot and 
dense matter created in nuclear collisions [10] . Therefore 
these photons were proposed to verify the existence of 
the QGP phase [TTJ [T^] . Spectra of thermal photons de- 
pend upon the fireball temperature and they can be cal- 
culated from the scattering cross-section of the processes 
like qq — ¥ gj, bremsstrahlung etc. Time evolution of the 
temperature can be calculated using hydrodynamics with 
appropriate initial conditions. Thus the spectra depend 
upon the equation of state (EoS) of the medium and they 
may be useful in finding a signature of the quark-gluon 
plasma [T2HH]. Recently thermal photons are proposed 
as a tool to measure the shear viscosity of the strongly 
interacting matter produced in the collisions [TTJ [TB]. 

Understanding the shear viscosity of QGP is one of the 
most intriguing aspects of the experiments at Relativistic 
Heavy Ion Collider (RHIC). Analysis of the experimental 
data collected from RHIC shows that the strongly cou- 
pled matter produced in the collisions is not too much 
above the phase transition temperature T c and it may 
have extremely small value of shear viscosity r\. The ra- 
tio of the shear viscosity 77 to the entropy density s i.e. 
•q/s is around 1/4-7T which is the smallest for any known 
liquid in the nature [19]. In fact the arguments based on 
AdS-CFT suggest that the values of rj/s can not become 
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lower than l/Ait. This is now known as Kovtun-Son- 
Starinets or 'KSS- bound' [20] • Thus the quark-gluon 
plasma produced in RHIC experiments is believed to be 
in a form of the most perfect liquid [21]. No wonder ideal 
hydrodynamic appears to be the best description of such 
matter as suggested by comparison between the experi- 
mental data[22] and the calculations done using second- 
order relativistic hydrodynamics [23-30 . 

However there remain uncertainties in understanding 
the application and validity of the hydrodynamical proce- 
dure in relativistic heavy-ion collision experiments. It is 
only very recently realized that the effect of bulk viscos- 
ity can bring complications in the hydrodynamical de- 
scription of the heavy-ion collisions. Generally it was 
believed that the bulk viscosity, £ does not play a signif- 
icant role in the hydrodynamics of relativistic heavy-ion 
collisions. It was argued that that since C scaled like 
e — 3P at very high energy the bulk viscosity may not 
play a significant role because the matter might be follow- 
ing the ideal gas type equation of state [31]. But during 
its course of expansion the fireball temperature can ap- 
proach values close to T c . Recent lattice QCD results 
show that the quark-gluon matter do not satisfy ideal 
EoS near T c and the ratio Q/s show a strong peak around 
T c [32] [33]. The bulk viscosity contribution in this regime 
can be much larger than that of the the shear viscosity. 
Recently the role of bulk viscosity in heating and expan- 
sion of the fireball was analyzed using one dimensional 
hydrodynamics [33]. Another complication that bulk vis- 
cosity brings in hydrodynamics of heavy-ion collisions is 
phenomenon of cavitation|35j. Cavitation arises when 
the fluid pressure becomes smaller than the vapour pres- 
sure. Since the bulk viscosity (and also shear viscosity) 
contributes to the pressure gradient with a negative con- 
tribution, it may be possible for the effective fluid pres- 
sure to become zero. Once the cavitation sets in, the 
hydrodynamical description breaks down. It was shown 
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in Rcf. [35) that cavitation may happen in RHIC exper- 
iments when the effect of bulk viscosity is included in 
manner consistent with the lattice results. It was shown 
that the cavitation may significantly reduce the time of 
hydrodynamical evolution. 

Keeping the above discussion in mind, we aim to study 
the effect of bulk viscosity and cavitation on the thermal 
photon production in heavy ion collisions. As far as we 
know no such study exists in the literature. Further- 
more, the calculations for the photon production rates 
in the absence of viscous effects, are done using thermal 
distribution function of the particle species (e.g., quark, 
anti-quark etc.)[S]. However, it is well known that vis- 
cous effects can lead to the modification of the thermal 
distribution functions [3B]. This may have observational 
effect on the photon spectraJTS] . In this work we incorpo- 
rate the viscous modification in the distribution function 
arising due to bulk and shear viscosities. Finite ( effect 
can either significantly reduce the time for the hydro- 
dynamical evolution (by onset of cavitation) or it can 
increase the time by which the system reaches T c . More- 
over the non-ideal gas EoS can also significantly influence 
the hydrodynamics. In what follows, we use equations of 
relativistic second order hydrodynamics to incorporate 
the effects of finite viscosity. We take the value of £/s 
same as that in Ref. [35] and keep rj/s = l/4ir. Further 
we use one dimensional boost invariant hydrodynamics 
in the same spirit as in Refs.[3H 131]. One of the limi- 
tations of this approach is that the effects of transverse 
flow cannot be incorporated. As the boost-invariant hy- 
drodynamics is known to lead to under-estimation of the 
effects of bulk viscosity [33], we believe that our study of 
the photon spectra will provide a conservative estimate 
of the effect. However it should also be noted that the 
effect of transverse flow could remain small as cavitation 
can restrict the time for hydrodynamical evolution. 



II. FORMALISM 



A. Viscous Hydrodynamics 



contribution. The corresponding hydrodynamics equa- 
tions are given by, 



De + (e 
(e + P) Du a - 



-P)0 
V a P 



(3) 
(4) 



where D 
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V Q = Auad^ and 



A v gives the symmetrization. 



We employ Bjorken's prescription [37] to describe the 
one dimensional boost invariant expanding flow, were we 
use the convenient parametrization of the coordinates us- 
ing the proper time r = y/t 2 — z 2 and space-time rapidity 
y = h ^[fzfji t — t cosh y and z = r sinhy. Then the 
four velocity is given by, 

it' 1 = (cosh y , 0, 0, sinh y) . (5) 

We note that with this transformation of the coordinates, 
D = #- and 6 = 1/t. 

Form of the energy momentum tensor in the local rest 
frame of the fireball is then given by |38H41] : 



( e 

Pj_ 

P ± 

\ P z 



(6) 



where the effective pressure of the expanding fluid in the 
transverse and longitudinal directions are respectively 
given by 



Pz=P 



P + U+ 



(7) 



Here <E> and n are the non-equilibrium contributions to 
the equilibrium pressure P coming from shear and bulk 
viscosities respectively. Respecting the symmetries in the 
transverse directions the traceless shear tensor has the 
form = diag($/2, $/2, -$). 

In the first order Navier-Stokes dissipative hydrody- 
namics 



We represent the energy momentum tensor of the dis- 
sipative QGP formed in high energy nuclear collisions as 



PA' 



where e, P and w' 1 are the energy density, pressure and 
four velocity of the fluid element respectively. The oper- 
ator A' 1 " = — u v acts as a projection perpendic- 
ular to four velocity. The viscous contributions to 
are represented by 



nv 



AM"n 



(2) 



where 7r' tl/ , the traceless part of n M "; gives the contri- 
bution of shear viscosity and n gives the bulk viscosity 



n 



-(d^u^ and 



nv 



(8) 



with £, n > and V^it,,) = 2V( M u„) — | A MZ ,V Q u Q . So 



^ for first order theories with Bjorken flow we have 



11=-- and . 
r 3r 



(9) 



The Navier-Stokes hydrodynamics is known to have in- 
stabilities and acausal behaviours [?2] [53]- second order 
theories removes such unphysical artifacts. 

We use causal dissipative second order hydrodynamics 
of Isreal-Strewart [44] to study the expanding plasma in 
the fireball. In this theory we have evolution equations 
for n and $ governed by their relaxation times Tn and 
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tv. We refer [15JIIS] for more details on the recent devel- 
opments in the theory and its application to rclativistic 
heavy ion collisions. 

Under these assumptions, the set of equations (i.e., 
equation of motion (|3| and relaxation equations for vis- 
cous terms) dictating the longitudinal expansion of the 
medium are given by (40] |43j [47] 
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(11) 
(12) 



where $ = tt 00 — tt zz . The terms in the square bracket 
in Equation( |ll[ ) are needed for the conformality of the 
theory[48 . The coefficients (3q and /?2 are related with 
the relaxation time by 



Tn = C A) ) Tir = 2 *?#2< 



(13) 



We use the Af = 4 supersymmetric Yang-Mills theory 
expressions for tv and Ai [48^50] : 



and 



Ai 



2 - In 2 
2vrT 



2ttT 



(14) 



(15) 



We set r OT (T) = Tn(T) as we don't have any reliable pre- 
diction for tt-t[3~4"]. 

In order to close the hydrodynamical evolution equa- 
tions (10 - 12 1 we need to supply the EoS. 



B. Equation of state, f/s and rj/s 

We are interested in the effect of bulk viscosity on the 
hydrodynamical evolution of the plasma and recent stud- 
ies show that near the critical temperature T c effect of 
bulk viscosity becomes important [5 ll - f55] . We use the re- 
cent lattice QCD result of A. Bazavov et aZ.[35] for equi- 
librium equation of state (EoS) (non-ideal: e — 3P ^= 0). 
Parametrised form of their result for trace anomaly is 
given by 



3P 



J>4 



V 



1 + cxp 



y2 



(I4 



(16) 

where values of the coefficients are d,2 = 0.24 GeV 2 , dn = 
0.0054 GeV 4 , c x = 0.2073 GeV, and c 2 = 0.0172 GeV [35]. 
Their calculations predict a crossover from QGP to 
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FIG. 1. (e - 3P)/T 4 ,C/s (and r//s = l/4vr) as functions 
of temperature T. One can see around critical temperature 
(T c = .190 GeV) C 3> f) and departure of equation of state 
from ideal case is large. 



hadron gas around 0.2-0.18 GeV. We take critical tem- 
perature T c as 0.19 GeV throughout the analysis. The 
functional form of the pressure is given by [32] 



P(T) P(T ) 
T 4 T* 



dT' 



T 



e-3P 
y5 



(17) 



with T = 50 MeV and P(T ) = [35J. 

From Equations ( 16 ) and p7] ) we get e and P in terms 
ofT. 

We rely upon the lattice QCD calculation results for 
determining £/s. We use the result of Meyer [33 , which 
indicate the existence a peak of Q/s near T c , however 
the height and width of this curve are not well under- 
stood. We follow parametrization of Meyer's result from 
Ref . [55] . given by 



aexp 



T c -T 
AT 



+ b 



for T > T c , (18) 



where b = 0.061. The parameter a controls the height and 
AT controls the width of the Q/ s curve and are given by 



0.901, AT 



14.5' 



(19) 



We will change these values to explore the various cases 
of (/s to account for the uncertainty of the height and 
width of the curve. 

We use the lower bound of the shear viscosity to en- 
tropy density ratio known as KSS bound 20J 



n/s = 1/4tt 



(20) 



in our calculations. We note that the entropy density is 
obtained from the relation 
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FIG. 2. Various bulk viscosity scenarios by changing the 
width of the curve through the parameter AT. 



renchc et at showed that two loop level bremsstrahlung 
process contribution to photon production is as impor- 
tant as Compton or annihilation contributions evaluated 
up to one loop level [55], They also discussed a new mech- 
anism for hard photon production through the annihila- 
tion of an off-mass shell quark and an antiquark, with 
the off-mass shell quark coming from scattering with an- 
other quark or gluon. These processes in the context 
of hydrodynamics of heavy ion collisions were studied 
in Refs.[13j [14]. Until recently only the processes of 
Compton scattering and qq- annihilation were considered 
in studying the photon production rates. 

The production rate for hard (E > T) thermal pho- 
tons from equilibrated QGP evaluated to the one loop 
order using perturbartive thermal QCD based on hard 
thermal loop (HTL) resummation to account medium ef- 
fects. The Compton scattering and qq- annihilation con- 
tribution to the photon production rate is[IJ[2l[5] 



In Fig. T] we plot the trace anomaly (e — 3P) /T 4 and 
C/s for desired temperature range. We also plot the 
constant value of n/s = 1/A-w for a comparison. It is 
clear that the non-ideal EoS deviates from the ideal case 
(e = 3P) significantly around the critical temperature. 
Around same temperature £/s starts to dominate over 
n/s significantly. We would like to note that these re- 
sults are qualitatively in agreement with Ref . |34j . 

In Fig. [2] we show the change in bulk viscosity profile 
by varying the width of the Q/s curve by keeping the 
height intact. 
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cE N 
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(22) 

where the constant c rs 0.23 and a and a s are the elec- 
tromagnetic and strong coupling constants respectively. 
In the summation / is over the flavours of the quarks 
and e/ is the electric charge of the quark in units of the 
charge of the electron. 

The rate of photon production due to Bremsstrahlung 
processes is given by [56 



C. Thermal photons 

During QGP phase thermal photons are originated 
from various sources, like Compton scattering q{q)g — > 
q(q)"f and annihilation processes qq — > gj. Recently Au- 
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(23) 

where Jt ~ 1.11 and Jj, m 1.06 for two flavours and 
three colors of quarks[T4"]. The expression for I(E,T) is 
given by 



I(E,T) 



3C(3) 



7T 2 E (E 



6 T \T 



ln(2)+4Li 3 (-e-l E l/ T )+2^ 



E 



)li 2 (~e-l E l/ T )-(|) 2 m(l + e 



-|E|/T 



(24) 



where Li are the polylogarithmic functions given by scattering in the medium is given by, 

E _dN 

+°° „n d 

Lia(z) = £ 



n=l 



z" 
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d 4 xd 3 p 3n 5 
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£4] ETe~ E / T (J T -J L ). 

(25) 



We use the parametrization of a s (T) by Karsch[57]: 

()7T 



Next the rate due to qq- annihilation with an additional " s(T) = (33^2^) ln(8T/T c ) 
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FIG. 3. Hard thermal photon rates in QGP as a function of 
energy for a fixed temperature T=250 MeV. Photon rates are 
plotted for different relevant processes. 



for our rate calculations. Here Nt is the number of 
quark flavors in consideration. 

In Fig. [3], we plot the different photon rates for a 
fixed temperature T = 250 MeV. It shows the contri- 
butions from Bremsstrahlung (Brems), annihilation with 
scattering (A+S) and Compton scattering together with 
qq- annihilation (C+A). Bremsstrahlung contributes to 
the photon production rate upto E ~ I GeV only, 
afterwards A+S and C+A processes become dominant. 
We might mention here that this observation is in 
agreement with Ref. [14] . 

The total photon rate is obtained by adding different 
temperature depended photon rate expressions. Once 
the evolution of temperature is known from the hydro- 
dynamical model, the total photon spectrum is obtained 
by integrating the total rate over the space time history 
of the collision 1 58 , 



dN 
d 2 p T dy 



dN 



(27) 



= Q drr 

J T 



^ d 3 pd 4 x / 

, / dN 

where To and t\ are the initial and final values of time 
we are interested. y nU c is the rapidity of the nuclei 
whereas Q is its transverse cross-section. For a Au nu- 
cleus Q <~*j 180/m 2 . pt is the photon momentum in di- 
rection perpendicular to the collision axis. The quantity 

is Lorentz invariant and it is evaluated in the 



( TP dN \ 



local rest frame in equation ([27]). Now the photon en- 
ergy in this frame, i.e., in the frame comoving with the 
plasma, is given as pTCosh{y — y'). So once the rapidity 
and pt are given we get the total photon spectrum. 



D. Viscous corrections to the distribution functions 

Viscous effects contribute in two ways in kinetic theory: 
Firstly, it can change the width (temperature) of the dis- 
tribution function. Secondly, it can modify the momen- 
tum dependence of the distribution function. The first 
effect is incorporated when we calculate the temperature 
as a function of time using dissipative hydrodynamics. 
To include the second effect one needs to compute the 
change in the distribution function as a function of mo- 
mentum using the techniques of kinetic theory [36]. In 
the following we give some details of such a calculation. 

In section II C[ while writing the photon rates, we have 
used Boltzmann distribution function of type f = fa = 
e — pu/T _ j n orc [ er t incorporate the modification due 
to viscous effects we write the distribution function as 
/ = fo + Sf, with 5f = Sf v + Sf c , where Sf v and Sf c rep- 
resent change in the distribution function due to shear 
and bulk viscosity respectively. We calculate Sf using 
14-moment Grad's method. It ought to be mentioned 
that recent results show that calculation of Sf using this 
method fails near freezeout region by making Sf even 
larger than f and / < [SMrJT] . It is therfore important 
to here that we are applying these corrections to calcu- 
late the photon production rate of hard thermal photons 
in the regime T > T c . We have found that for pr be- 
low 3 GeV, this approximation is reasonable but beyond 
it, this approximation breaks down as the contribution 
arising from the viscous correction Sf to the distribution 
function becomes larger than / [55] ■ With this caveat we 
proceed to calculate Sf applying the techniques used in 
Refs.[12|3H]. 

We write the viscous correction to the (Boltzmann) 
distribution function as 

f(p) =fo + §f = fo + Sf n + Sf c (28) 

= f (i + v (a ^> + ^V^ej 

where we restrict the corrections to / upto quadratic or- 
der in momentum. In order to find the coefficients A and 
C we first express the energy momentum tensor using /, 
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Now from Eq. 30 we get the correction Sf v due to the 
shear viscosity as given in Ref. [35] by finding out C and 



G 



we will not repeat that calculation here. Next we will 
find out the coefficient A by constructing a fourth rank 
symmetric tensor out of A Miy and u M representing the 



term in square brackets in Eq. 31 
d 3 p 



A 
2T 3 



= a a (u^u v u a uP) (32) 



+Oi (A^ u u a u^ + permutations) 
+a 2 (A^A a0 + A ,ia A l/ ^ + A^A va ) 



Now substituting this expression in Eq. 31 and by not 
ing A M „ u v = , A^A"" = 3 and A""A MQ = A; 
get £ = 5a 2 . Now by contracting both sides of Eq. 
with i (A^A Q(3 + A^A^ + A^A^) we get, 
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Evaluating this expression in the local rest frame of the 
fluid = (1,0) we get 
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(2^)3 J B /o|PI 



(34) 



Now for a Boltzmann gas with fo = e~ pu / T we can cal- 
culate the integral and comparing the result with that of 
the entropy density s of an ideal boson gas[3H] we find, 
A = | C/s. So the viscous correction to the distribution 
function due to both shear j38] and bulk viscosities are 
given as 

/ = fo (l + ^VV^fl + | ^VA^ej . (35) 

Using the Bjorkcn's flow one can calculate V( Q u^) and 
A a pQ for the present problem. Four velocity can be writ- 
ten as u a = (cosh y', 0, 0, sink y'), where y' is the rapid- 
ity. Let the four momentum of a particle be parametrised 
as p a = (mTCOshy,pTCOs4> p ,pTsin(f)p,mTCOshy), where 
m T = v\ +w 2 [3"g]. 

Finally we write the distribution function including 
viscous correction as 



f = fo 1 



2T 3 
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5 2T 3 
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This expression for distribution function will be used 
in estimating the photon production rates in Eq. 27 . 



III. RESULTS AND DISCUSSION 

In order to understand the temporal evolution of 
temperature T(t), pressure P(t) and viscous stresses - 



TABLE I. Initial conditions for RHIC 



Vnuc T T() 

(fm/c) (GeV) 



5.3 0.5 



.310 



$(r) and LI(t), we numerically solve the hydrodynami- 
cal equations describing the longitudinal expansion of the 
plasma: Eqs. 10p2 . We use the 



-ideal EoS obtained 



from Eq. 16 and Eq. 17 



Information about viscosity 
coefficients £ and r] are obtained from Eqs. 18p0 using 
Eq. [2l] . We need to specify the initial conditions to solve 
the hydrodynamical equations, namely the initial time r 
and To. We use the initial values relevant for RHIC ex- 
periment given in Table I, taken from Ref.[13]. We will 
take initial values of viscous contributions as $(t ) = 
and II(to) = 0. We would like to note that our hydrody- 
namical results are in agreement with that of Ref . [35] . 

Once we get the temperature profile we calculate 
the photon production rates. Total photon spectrum 



E 



dN 



d3pd i x (as a function of rapidity, y and transverse mo- 
mentum of photon, pt) is obtained by adding different 

22|23|25 



photon rates using Eqs. 



and convoluting with 



the space time evolution of the heavy-ion collision with 
Eq. 27 . The final value of time T\ is the time at which 



temperature evolves to critical value 77 , i.e.; T{j\) = T c . 
In all calculations we will consider the photon production 
in mid-rapidity region (y = 0) only. 

We will be exploring various values of viscosity and its 
effect on the system. Since there is an ambiguity regard- 
ing the height and width of £/s curve, we will vary the 



parameters a and AT from its base value given in Eq. 19 



By this we will able to study the effect of variation of C 
on the system. The varied values of the parameters are 
represented by a' and AT'. We note that unless spec- 
ified we will be using the base values of bulk viscosity 
parameters (from Eq. 19 ) in our calculations. Through- 



out the analysis we will keep the shear viscosity 77 to its 



base value given by Eq. 20 



In order to understand the effect of non-ideal EoS in 
hydrodynamical evolution and subsequent photon spec- 
tra we compare these results with that of an ideal EoS 
(e = 3P). We consider the EoS of a rclativistic gas of 
massless quarks and gluons. The pressure of such a sys- 
tem is given by 



P = aT 4 



16 



21 
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90 



(37) 



where Nt = 2 in our calculations. Hydrodynamical evo- 
lution equations of such an EoS within ideal (without 
viscous effects) Bjorken flow can be solved analytically 
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FIG. 4. Temperature profile using massless (tdeal) and non- 
ideal EoS in RHIC scenario. Viscous effects are neglected 
in both cases. System evolving with non-ideal EoS takes a 
significantly larger time to reach T c as compared to ideal EoS 
scenario. 



and the temperature dependence is given by [57] 



T = Tn 



1/3 



(38) 



where tq and To are the initial time and temperature. 
While considering the viscous effect of this ideal EoS, 
we will solve the set of hydrodynamical equations flO] - 
111, since effect of bulk viscosity can be neglected in the 
relativistic limit when the equation of state P = e/3 is 
obeyed [31"] . 



Hydrodynamics with non-ideal and ideal EoS 

Fig. [i] shows plots of temperature versus time for the 
ideal and non-ideal equation of states. The temperature 
profiles are obtained from the hydrodynamics without 
incorporating the effect of viscosity. The figure shows 
system with non-ideal EoS takes almost the double time 
than the system with ideal massless EoS to reach T c . 
So even when the effect of viscosity is not considered, 
inclusion of the non-ideal EoS makes significant change 
in temperature profile of the system. This can affect the 
corresponding photon production rates. 

Next we analyse the viscous effects on the tempera- 
ture profile. The role of shear viscosity in the boost 
invariant hydrodynamics of heavy ion collisions, for a 
chemically nonequilibrated system, was already consid- 
ered in Ref.[T7]. We consider possible combinations of 
$ and II in non-ideal EoS case and study the corre- 
sponding temperature profiles as shown in Fig. IB]. As 
expected viscous effects is slowing down temperature evo- 
lution. For the case of non zero bulk and shear viscosities 



FIG. 5. Figure shows time evolution of temperature with non- 
ideal EoS for different combinations of bulk (II) and shear ($) 
viscosities. Non zero value of bulk viscosity refers to E qs.[18f 
1 1 91 and non zero shear viscosity is calculated from Eq. 1201. 




FIG. 6. Longitudinal pressure P z for various viscosity cases 
shown in Fig. [2]. 



(II ^ 0; $ 7^ 0), temperature takes the longest time to 
reach T c as indicated by the top most curve. This is 
about 35% larger than the case without viscosity (the 
lowest curve). The remaining two curves show that the 
bulk viscosity dominates over the shear viscosity when 
the value of T approaches T c and this makes the system 
to spend more time around T c . However the intersection 
point of the two curves may vary with values of a and 
AT as highlighted by Fig. [2]. 



Non-ideal EoS and Cavitation 

Let us note here the fact that, the bulk viscosity con- 
tribution LT is negative|35]. From the definition of Ion- 
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FIG. 7. Temperature is plotted as a function of time. With 
peak value (a) of (/s remains same while width (AT) varies. 
In all the three curves, solid lines end at cavitation time t c 
denoted by a dark circle. The dashed lines in each curves 
show how the system would evolve till T c if cavitation is ig- 
nored. Figure shows that larger the width parameter shorter 
the cavitation time. 



gitudinal pressure P z = P + II — $ it is clear that if 
either II or $ is large enough it can drive P z to nega- 
tive values. Pz = defines the condition for the onset 
of cavitation. During the course of expansion when P z 
vanishes, the fluid will break apart into fragments and 
the hydrodynamic treatment will become invalid (see for 
e.g., Ref. (35|). Recent experiments at RHIC suggest rj/s 
to its smallest value ~ 1/Att. Such a small value of ij/s 
alone is inadequate to induce cavitation. Therefore we 
vary the bulk viscosity values by changing a and AT to 
study the effect cavitation. In the discussion that follows 
we will use r c to denote the time when cavitation occurs. 

In Figs.[6][7] we plot P z and T as functions of the 
proper time for different values of AT while keeping a 
(=0.901) fixed. As may be inferred from Fig. [6], higher- 
value of AT leads to a shorter cavitation time. For the 



values of a and AT given by Eq. 19 we find that around 
t c = 2.5 fm/c, P z becomes zero as shown by the low- 
est curve in Fig. jfjj. In this case, the cavitation occurs 
when the temperature reaches the value about 210 MeV, 
as may be seen in Fig. [Fj. Had we ignored the cavitation, 
the system would have taken a time ry = 5.5 fm/c to 
reach T c , which is significantly larger than r c . This shows 
that cavitation occurs rather abruptly without giving any 
sign in the temperature profile of the system. The hydro- 
dynamic evolution without implementing the cavitation 
constraint can lead to over-estimation of the evolution 
time and the photon production. 

We have carried out a similar analysis shown in Figs. [6| 
by keeping AT fixed (= T c /14.5) and varying param- 
eter a. In Fig. [8] we show the cavitation times corre- 
sponding to changes in a and AT (denoted by a' and 



FIG. 8. Cavitation time r c as a function of different values of 
height (a') and width (AT') of £/s curve. 



AT'). The dashed curve in Fig. [8 shows t c as a function 
of a, while keeping AT fixed. The curve shows that r c 
decreases with with increasing a. Solid line shows how t c 
varies while keeping a fixed and changing AT. 



Thermal Photon Production 

We have already seen that the calculation of photon 
production rates require the initial time To, final time 
Ti and T(r). r\ and T(r) are determined from the hy- 
drodynamics. Generally t\ is taken as the time taken 
by the system to reach T c , i.e.; r/. Since hydrodynamics 
ceases to be valid beyond the cavitation time, we must set 
Ti = T c - Thus photon production from QGP will be in- 
fluenced by cavitation, temperature profile and non-ideal 
EoS near T c . 

Fig. [9] shows the photon production rate calculated us- 
ing ideal (massless) and non-ideal EoS. The figure shows 
that non-ideal EoS case can yield significantly larger pho- 
ton flux as compared to the ideal EoS. At px = 1 GeV, 
photon flux for the non-ideal EoS is about 60% larger 
than that of ideal EoS case. This is because the calcula- 
tion of the photon flux is done by performing time inte- 
gral over the interval between the initial time To and the 
final time ti . To is same for both the system while the Ti 
for the case with non-ideal EoS is two times larger than 
the ideal EoS, as may be seen in Fig. El. Further, since 
the non-ideal EoS allows the system to have consistently 
higher temperature over a longer period as compared to 
the massless ideal-gas EoS, more photons are produced. 

Next, we consider the question that how the cavita- 
tion can affect the photon production rate. We empha- 
size that the rates should only be integrated up to the 
cavitation time t c . Fig. 10 shows the case when there is 
no viscous correction to the distribution function. In the 
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FIG. 9. Photon flux as function of transverse momentum for 
different equation of states. No effect of viscosity included in 
the hydrodynamical equations and in the distribution func- 
tions. 




FIG. 10. Photon spectrum obtained by considering the effect 
of cavitation (dashed line). For a comparison we plot the 
spectrum without incorporating the effect of cavitation (solid 
line). Viscous correction to the distribution function is not 
considered. 



dashed curve the effect of cavitation is taken into account 
and T\ = t c = 2.5 fm/c. The solid line represents the 
same case but without the effect of the cavitation and 
Ti = tj — 5.5 fm/c. It can be seen from the curve that 
ignoring cavitation leads to an over-estimation of the rate 
by about 200% at p T = 0.5 GeV and about 50% at p T = 2 
GeV. It is thus clear that the information about the cav- 
itation time is crucial for correctly estimating thermal 
photon production rate. 

Fig. 11 shows the similar comparison between cavita- 
tion and no-cavitation cases as in Fig. 10 , but with the 
inclusion of viscous correction to the distribution func- 
tion. The solid curve shows the case when cavitation 



FIG. 11. Same as in Fig.[lO] but with incorporating viscous 
correction to the distribution function. 



is ignored. The dashed curve shows the effect of cavi- 
tation. First we note that as can be seen from Fig. [I], 
C/s 3> rj/s near T c , while rj/s dominates over C,/s when 
T T c . If one ignores the cavitation effect, then the hy- 
drodynamical code can allow for temperature to evolve 
upto T c . Moreover, we have also observed that 8fe con- 
tribution dominates over 5f v contribution for px <1.5 
GeV (A similar behaviour is reported in the Ref. |62|). 
Let us first note that for Sf = 0, the photon flux without 
cavitation is higher as compared the same with cavita- 
tion at the starting p T (~ 0.5 GeV). This feature also 

The 



continues when Sf ^ as can be seen in Fig. 11 
negative contribution of 5fe on the curve without cavi- 
tation (solid curve) makes it plummate faster than the 
curve with cavitation with increasing px (dashed curve) 
and both the curves intersect at pr ~1.1 GeV. Beyond 
this, the shear viscosity starts becoming more effective 
and prevents this faster plummation of the solid curve. 
Moreover, 6fc contribution dominates at pt < 1.5 GeV 
Thus in the regime pt < 1.5, ignoring the cavitation ef- 
fect can lead to over-estimation of the photon flux; e.g., at 
Pt = 0.5 GeV over-estimation is 150%. Whereas there is 
an under-estimation of the photon flux when the effect of 
cavitation is not included in the higher regime. E.g., for 
the following pt values 1.5, 2.0 and 3.0 (in GeV) under- 
estimations are around 50%, 65% and 30% respectively. 

In Fig. |12] we plot photon production rates for various 
cavitation times obtained by varying AT (with a = 0.901 
is fixed) . Here the enhancement in the photon production 
when AT is reduced to half of its base value is about 75% 
at p T = 0.5 GeV and about 55% at p T = 1 GeV. A fur- 
ther reduction of the parameter value to AT/4 is enhanc- 
ing the photon production by about 120% at pt = 0.5 
GeV and about 85% at pt = 1 GeV. The reason is a re- 
duction in AT amounts to increase in the cavitation time 
(see e.g., Fig. [6]), which in turn increases the time inter- 
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FIG. 12. Photon production rates showing the effect of dif- 
ferent cavitation time. Viscous corrections to the distribution 
functions has been included. 



and a (=0.901). The solid curve shows the case Sf = 0. 
We can see from the figure that when we include only 
the shear viscosity correction (Sf = Sf v ), it enhances 
the photon production particularly in high pj< regime [IB]. 
Next we consider the corrections arising due to bulk vis- 
cosity Sf = Sfc; only. The effective pressure due to II is 
negative and we have seen that it is crucial for the cavita- 
tion to occur (See Eq. [7] ) . As a result number of particle 
with higher momenta are decreasing and therfore there is 
a reduction in photon rate as compared to the case with- 
out the viscous correction [62], As it is clear from the 
figure, the effect of bulk viscosity is to oppose the contri- 
bution from the shear viscosity. The combined effect of 
both shear and bulk viscosity corrections (Sf = Sf^+Sf^) 
can be seen in the graph just below the top-most curve. 
It is clear from the graph that effect of viscous correc- 
tions to the distribution functions increases in the high 
Pt regime. 
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FIG. 13. Viscous corrections to the distribution function and 
photon production rate. The solid line shows the photon pro- 
duction rate without the viscous corrections to the distribu- 
tion function and the other lines shows the cases with addition 
of viscous corrections due to shear and bulk viscosities. 



val over which photon production is calculated. There- 
fore this increases the photon flux. 

In Fig. [13] we show the effect of viscous corrections 
to the distribution function on photon production. Here 
we consider the cavitation scenario with AT (= T c /14.5) 



IV. SUMMARY AND CONCLUSIONS 

Using the second order relativistic hydrodynamics we 
have analyzed the role of non-ideal effects near T c arising 
due to the equation of state, bulk-viscosity and cavita- 
tion on the thermal photon production from QGP. Since 
the experiments at RHIC imply extremely small values 
for 77/s, we take the value 1/47T for shear viscosity. We 
have shown using non-ideal EoS from the recent lattice 
results that the hydrodynamical expansion gets signifi- 
cantly slow down as compared to the case with the mass- 
less EoS. This, in turn, enhances the flux of hard thermal 
photons. 

Bulk viscosity plays a dual role in heavy-ion collisions: 
On one hand it enhances the time by which the system 
attains the critical temperature, while on the other hand 
it can make the hydrodynamical treatment invalid much 
before it reaches T c . Another result we would like to em- 
phasize is that if the viscous correction Sf to the distri- 
bution function is not included, then ignoring cavitation 
can lead to a significant over-estimation of the photon 
production rate. But when the viscous corrections are 
included the situation can become more complex. In the 
low px regime (<1.1 GeV) there is a significant over- 
estimation in the photon flux if cavitation is ignored. On 
the other hand, in the high pt regime (>1.1 GeV) the 
photon flux is under- estimated when the cavitation is ig- 
nored! 
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